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Abstract 

Primordial magnetic fields of cosmologically interesting field strengths can be generated from grav- 
itationally coupled electrodynamics during inflation. As the cosmological constraints require this 
to be power law inflation it is not possible to generate at the same time the curvature perturbation 
from inflation. Therefore here a completion is considered whereby the large scale magnetic field is 
generated during inflation and the primordial curvature mode in a subsequent era from a curvaton 
field. It is found that constraints on the model to obtain strong magnetic fields and those to sup- 
press the amplitude of the curvature perturbation generated during inflation can be simultaneously 
satisfied for mag netic seed fields B s ~ 10" 30 G. 



1 Introduction 

Observations indicate the existence of large scale magnetic fields in the universe. Evidence for 
magnetic fields at galactic and cluster scale has been obtained for decades with field strengths of 
order 10 -6 G [I]. However, over recent years several groups have found indications of truly cosmo- 
logically magnetic fields, pervading space with intensities of order 10 -15 G [2]. There is a multitude 
of proposals to explain the existence of cosmological magnetic fields (for recent reviews, see, e.g., 
[3]). Placing the generation mechanism in the very early universe during inflation by amplifying 
significantly perturbations in the electromagnetic field requires to break conformal invariance in 
flat backgrounds [4] which however is not the case in open geometries [5] (see, however, [6]). 

The conformal invariance of Maxwell's equations in four dimensions is broken in models in 
which the electromagnetic field is gravitationally coupled. Couplings between curvature terms 
and the Maxwell tensor in the Lagrangian of the form RF flll F llI/ , R VLV F ilR F v R and R iJLU \ K F tll/ F Xk 
are present when describing, e.g., the propagation of a photon in a curved background [7J. In 
Fourier space the resulting mode equation for the gauge potential of the electromagnetic field in an 
expanding background includes a term which during inflation can cause sufficient amplification on 
superhorizon scales so as to generate a strong enough magnetic field in order to seed the galactic 
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dynamo. This was first proposed in [5] and investigated in more detail in |S]. 

In [8] the final spectrum of the magnetic field is calculated assuming the inflationary stage is 
directly matched to the standard radiation dominated era at some conformal time t = t\, so that 
the scale factor has the form 



a(r) 



and the line element has the form ds 2 




(1.1) 



(t)(—cIt 2 + dx 2 + dy 2 + dz 2 ). In the following a\ = 1. 



For t < ti, de Sitter inflation takes place for j3 = — 1 and power law inflation for 
In [8] the fractional magnetic field energy density, r 
been found as 



-oo < j3 < -1. 
^ on a galactic scale of order 1 Mpc has 



r(u G ) = 10 



-79+52, [r{ly)] 2 ( 1 
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where z/ = |/3 + ||, £2 = yw~^o an d -f^i is the Hubble parameter at the beginning of the radiation 
dominated era, at r = t\ which determines the reheat temperature. It was found that a magnetic 
field strength B > 1CP 20 G corresponding to r > 1CP 37 can be reached for (3 < —2.8 and < 
10~ 18 , depending on j3 (cf. figure [Q). The constraint on /3 prevents that the observed curvature 
perturbation is generated during inflation since the resulting spectral index is too small. In slow 
roll inflation power law inflation [S] is realized by an exponential potential of the inflaton <p of the 
form (e.g, [TO] ) 



V{ct>) = V % exp 



p M P 



(1.3) 



where the index i indicates initial values and p determines the scale factor in cosmic time dt = adr, 
that is a ~ t p . Thus p = -J^- The slow roll parameters are given by e = ^ = -. Ending inflation by 
bubble formation puts the constraint p < 10 |11| 112]. p as a function of f3 is monotonically growing 
for j3 < 0, approaching 1 for /3 — > —00 and at j3 = —2.8, p = 1.56, so that the upper limit on p 
does not put any additional constraint on the model at hand. The spectral index of the curvature 
perturbation created during inflation, n s = 1 + 2r/ — 6e, is in the range —1 and —0.29 for f3 < —2.8. 
The power spectrum of the curvature perturbation evaluated at horizon crossing is given by [10] 
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where /i = | + . The amplitude of the curvature perturbation generated during inflation should 
only contribute negligibly to the final curvature perturbation. Starting with the slow roll equations 



3H<i> + 

H 2 



the evolution of the inflaton d> is found to be 
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Thus the potential is given by 
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This allows to calculate Hk = H((j)k) = H\k= a H- Using that aH = & and a\ = a{r\) = 1 it follows 
that 
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So that finally, using 
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where /c p is going to be chosen to be the pivot wavenumber of WMAP 7 today, k p = 0.002 Mpc -1 
[T3] . Thus the amplitude of the curvature spectrum at k p is determined by 
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In figure [T] a contour plot of log 10 'P ( f is shown. As can be appreciated from figure Q] satisfying the 
constraints to generate a magnetic field with B s ~ 10 -20 G corresponding to r ~ 10 -37 , that is 
f3 < —2.8 leads to curvature perturbations which are in general, too large, not only for the curvaton 
mechansim to work, but even larger than the observed amplitude = 2 A3 x 10 -9 [13]. However, 
it is possible to generate a magnetic field with B s ~ 10 -30 G corresponding to r ~ 10~ 57 which is 
the lower limit required to seed the galactic dynamo when the non vanishing cosmological constant 
is taken into account |14j . In order for the curvature perturbation created during inflation to be 
subdominant /3 and jjy 3 - have to be in the range, —4.8 < (3 < —3.6 and 10 -40 < < 10~ 34 which 
corresponds to reheat temperatures between 0.1 and 100 GeV. 



2 Generating the curvature perturbation 

In the curvaton model the curvature perturbation is generated after inflation has ended, by the 
curvaton field a which has been subdominant during inflation [15] . The spectral index of the final 
curvature perturbation is given by |llj 

7^ = 1 + 27^-26, (2.1) 

where 

Mp 8 2 V . . 

= a* (2 - 2) 
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Figure 1: Contour plot of log 10 7^ shown for k p = 0.002 Mpc . The contour lines corresponding 

to vf = 10~ 9 (thick, black dashed line) and vf = 10~ 12 (thick, blue long dashed line) are marked. 
The red, dot-dashed and green, dotted lines show the maximal value of the logarithm of the ratio 
of magnetic to background radiation energy density of the model [Sj, corresponding to r = 10~ 37 
(red, dot-dashed) and r = 10 -57 (green, dotted). 



and Mp is the reduced Planck mass Mp / (8tv) . Assuming the simplest model |X5|, [T6] the potential 
during inflation is defined to be 



V(<f>, a) = Vi exp 



p M P 
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where the contribution from the curvaton has to be sub leading to the one of the inflaton. Equation 
(|2.ip requires rjo-o- to be of the order of e for a nearly scale invariant spectrum. Thus 
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calculated at the time of first horizon crossing during inflation of the mode k = aH. So that the 
mass of the curvaton is determined by 
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This is shown in figure [2] for the region allowed by the constraints that Vf < 10 and r > 10 
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Figure 2: Contour plot of log 10 (j^j shown for k p = 0.002 Mpc 1 in the region allowed by the 
constraints < W~ 12 and r > 10~ 57 for = 2 A3 x 10 -9 and 1 - n s = 0.037 [13]. 
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The evolution of the background curvaton is determined by 

a + 3Ha + V a = 0. 



(2.6) 



It is assumed [TTJ that \V aa \ -C if 2 which in the model at hand results in m 2 t 2 <C 1 since p is 0(1) 
which implies that during and after inflation, before the curvaton becomes effectively massive, a 
stays approximately constant, a ~ er*. 

In the subsequent analysis we follow [16 . The end of inflation takes place at r = n when the 
Hubble parameter has the value H±. During the following radiation dominated era it evolves as 
H 2 ~ 1/a 4 so that 



The curvaton becomes massive at m 2 = H 2 so that 



(2.7) 
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To prevent an additional stage of inflation driven by the curvaton thus requiring the universe to be 
radiation dominated at T mass imposes p r ad{ T mass) ^™^») where a* is the value of the curvaton 
during inflation when the modes of the observable perturbations leave the horizon, assuming that, 
during the late stages of slow roll inflation, the change in a can be neglected. Then together with 
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the value of p ra d at the beginning of the radiation dominated era, p ra d{ T i) = 3MpH 2 , it follows 
that 

« -ur (2 ' 9) 

which is the same constraint as in the chaotic inflation model of |16j . 

In [15] two separate cases have been considered. If the curvaton decays during the radiation 
dominated era, that is until its decay it stays subdominant, the resulting curvature perturbation is 
given by [15] 

V C = ^f4\ ( 2 - 10 ) 

16 TT (J* 

using that cr& = <?* is approximately constant during inflation and S = -^ s —. Moreover, here 
Pa = ^m 2 a 2 . As in [16] it is found that at the time of decay, 

2 

Q _ a * a decay O -\-\\ 

Sdecay- ^ , (-IM 



1 p u mass 



and defining the decay constant of the curvaton T a , together with T a = Hdecay results in a<fe ca2/ /a r , 
(m a /Ta.) 1 / 2 implying 



9 / \ - 
at ( TTln > 2 



Sdecay ~6M p \ rj • (2,12) 

So that finally, the spectrum of the curvature perturbation, for S 'decay < 1, is found to be 



p ^^m;J [w p ){w p ) [w P ) ■ (2 - 13) 

Moreover, for the perturbations to be Gaussian H^/a* <C 1 [15]. This together with S 'decay < 1 
determines the range of possible values of a^/H^ to be 

k |r < -ir. (2.i4) 

Finally, the decay constant IV is determined by 

f^ = _L(M(£LVW'. (2 . 15) 

VAf Py / 9V ( \MpJ \H k ) \M p ) 1 ' 

The decay constant T a and the value of cr* is shown for different values of jf^ in figure [3] As 
can be seen T a /Mp = Hdecay > 10~ 40 which is the minimal value to ensure standard primordial 
nucleosynthesis . 

In the opposite case, Sdecay > 1> the curvaton dominates before decay and [15] 
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log 10 (IV/Mp), crJH k = VlO" 
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Figure 3: Upper panel: Contour plot of log 10 i^f~^- Lower panel: Contour plot of the corresponding 
values of log 10 ^jf~^j ■ All plots are shown for k p = 0.002 Mpc -1 in the region allowed by the 

constraints < 10 -12 and r > 10~ 57 for different values of (jj^V The curvature perturbation is 

assumed to be given by the best fit parameters of WMAP 7, V ( = 2.43 x 10~ 9 and 1 - n s = 0.037 
. It is assumed that the curvaton decays during the radiation dominated era, so that S 'decay < 1 • 
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In this case the Gaussianity requirement is always satisfied if the curvature perturbation are of the 
observed magnitude. Moreover cr* is completely determined by equation (|2.16p . 
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which is shown in figured] for the best fit values of WMAP 7 [13]. In this case S 'decay > 1 which 
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Figure 4: Contour plot of log 10 (tOp) shown for fc p = 0.002 Mpc 1 in the region allowed by the 
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that the curvaton dominates before decay, so that Sd ecay > 1. 

results in the constraint, |16j 
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This implies (mp) ^ ^ 4 ° wmcn ^ s satisfied in the allowed region in parameter space as 

can be seen from figures [2] and 01 



3 Conclusions 

In [5] it was shown that cosmologically relevant magnetic fields can result from gravitationally 
coupled electrodynamics motivated by the form of the one-loop effective action of the vacuum 
polarization in QED in a gravitational background. Perturbations in the electromagnetic field are 
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amplified during power law inflation determined by the exponent j3 and the value of the Hubble 
parameter at the end of inflation H\ which here also determines the reheat temperature. The 
ratio r of magnetic field energy density over background radiation energy density has to be larger 
than 10 -37 , corresponding to B s ~ 10 -20 G to seed the galactic dynamo in a universe with no 
cosmological constant. This value is reduced to r ~ 10 -57 and B s ~ 10 -30 G in a universe with 



In [8] it was found that there is a region in the (/?, y jj^J )-plane for which seed magnetic field 

are obtained with B s > 10 -20 G corresponding to r > 1CP 37 . However, the curvature perturbations 
generated during inflation for those values are incompatible with the observed nearly scale invariant 
spectrum. Therefore, here the possibility of completing this model of magnetic field generation with 
the curvaton mechanism has been considered. Thereby generating the curvature perturbations after 
inflation. Hence imposing that the amplitude of the curvature perturbations during inflation is 
negligible. This restricts significantly the parameter space in (3 and (^jjr^J allowing the creation of 

primordial magnetic fields only with amplitude B s ~ 10 -30 G, however, still satisfying the weaker 
constraint r > 1CP 57 . 

We have assumed that power law inflation is realized within slow roll inflation with an expo- 
nential potential. The curvaton is described by a simple quadratic potential. Assuming that the 
curvature perturbation due to the curvaton is determined by the best fit values of WMAP 7 and 
that the contribution of the curvature perturbation due to the inflaton is less than 10~ 3 of that 
of the curvaton the parameter space of the curvaton model has been explored. The corresponding 

values in the ^-|^-J)-plane have been found for the mass of the curvaton m a , its field value 
during inflation <r* and the decay constant in the two cases where the curvaton decays during 
radiation domination or during curvaton domination. 
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